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This  report  contains  a theoretical  study  of  the  electrical  corona  on  a long  wire. 
The  physics  of  the  phenomenon  Is  discussed  In  detail,  and  formulated  mathemati- 
cally In  terms  of  coupled  nonlinear  partial  differential  equations.  A number  of 
boundary- value  problems  for  these  equations  are  solved  to  predict  the  corona 
breakdown  field  and  the  d-c  corona  voltage-current  characteristic.  The  results 
are  In  good  agreement  with  experiment.  Finally,  the  d-c  corona  effect  on  large 
signal  propagation  along  a single  wire  above  ground  Is  examined  by  solving  a 
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lonlinear  transmission  line  equation.  The  purpose  of  this  study  is  to  identify 
the  set  of  basic  corona  equations  as  well  as  the  mathematical  nature  of  the  prob- 
lem. The  application  of  these  equations  to  analyze  the  corona  effect  on  aircraft 
long-wire  antennas  will  entail  the  solution  of  certain  time-dependent  problems  of 
:onsiderable  mathematical  complexity.  These  latter  problems  are  not  covered  in 
the  present  study,  and  should  form  the  subject  of  future  investigations. 
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SECTION  1 


INTRODUCTION 


On  the  surface  of  a good  conductor  the  electric  charges  are  highly 
concentrated  In  regions  of  small  radii  of  curvature,  such  as  tips,  wires, 
edges  and  comers.  As  a consequence  the  local  electric  field  can  attain 
an  enormous  Intensity.  It  can  easily  exceed  a certain  critical  value 
beyond  which  Ionization  or  breakdown  of  the  surrounding  air  sets  In.  This 
critical  or  breakdown  field  strength  Is  strongly  dependent  on  atmospheric 
condition,  and  on  conductor  geometry  and  dimensions.  At  sea  level  Its  value 
Is  typically  of  the  order  of  30  kilovolts  per  centimeter,  which  Is  by  no 
means  altogether  too  high  for  common  occurrence  of  the  breakdown  phenomenon. 
When  Ionization  Is  Initiated,  the  air  becomes  conducting  and  a discharging 
current  flows  through  It  from  the  conductor.  The  Ionization  and  the  dis- 
charge are  accompanied  by  a visual  display  In  the  form  of  an  eerie  bluish 
or  reddish  glow  crowning  the  hlgh-fleld  region.  The  discharge  has  thus 
acquired  the  name  of  a corona.  When  the  corona  appeared  on  top  of  the 
masts  of  a ship  In  a storm.  It  was  often  referred  to  by  the  ancient  mariners 
as  St.  Elmo's  fire.  The  corona  Is  maintained  until  the  discharge  brings  the 
local  field  strength  down  to  below  the  critical  value. 

The  corona  discharge  turns  out  to  be  an  extraordinarily  complex  phenomenon 
In  which  scores  of  processes  can  participate.  The  net  effect,  however.  Is 
Invariably  the  dissipation  of  electrical  energy  Into  heat,  light,  sound,  and 
chemical  energy.  Not  surprisingly,  early  study  of  the  corona  was  carried  out 
by  the  electrical  power  engineers.  On  a high-voltage  power  transmission  line 
the  surface  electric  Intensity  easily  exceeds  the  corona  onset  value.  The 
ensuing  corona  discharge  constitutes  a leakage  current  and  a power  loss.  For 
a power  line  going  over  long  distances  the  accumulative  corona  loss  can  be 
considerable.  Out  of  this  concern  a number  of  Investigators,  notably  Peek,^ 
have  conducted  a series  of  detailed  measurements  on  the  corona  properties,  such 
as  the  breakdown  field  rtrength,  the  power  loss,  and  the  dependence  on  atmos- 
pheric and  conductor  conditions.  The  result  Is  a set  of  empirical  formulas 
directly  applicable  to  power  transmission  calculations. 
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Another  area  of  electrical  engineering  In  which  the  corona  discharge 
can  exert  considerable  Influence  is  the  operation  of  large  linear  antennas. 

On  Air  Force  aircraft,  examples  of  such  antennas  are  the  LF/VLF  tralllng- 
wire  antenna  and  the  HF  fixed-wire  antenna.  The  transmitters  can  con- 
ceivably drive  these  antennas  Into  the  corona  region,  not  to  mention  such 
Intense  electric  field  sources  as  lightning  and  the  nuclear  electromagnetic 
pulse.  The  appearance  of  coronas  results  in  power  dissipation,  signal  dis- 
tortion and  noise  Interference.  Unfortunately  there  exist  as  yet  few  or 
no  empirical  formulas  for  a quantitative  estimation  of  these  effects,  un- 
like In  the  case  of  power  transmission  engineering.  This  is  due  to  the 
difficulty  of  making  direct  corona  measurements  on  the  antennas  of  a high- 
flying aircraft.  On  the  other  hand  the  hlgh-altltude  atmospheric  environment 
ind  the  free-flight  aircraft-antenna  configuration  are  not  easy  to  simulate 
or  measurement  on  the  ground.  It  Is  desirable  to  conduct  a theoretical 
dy  of  the  corona  discharge  from  basic  physical  principles. 

The  corona  discharge  Is  a particular  mode  of  the  electric  gas  discharge 
which  has  been  investigated  extensively  since  the  latter  half  of  the  19th 
century.  Out  of  this  long-standing  effort  came  the  discovery  of  the  electron, 
the  invention  of  the  X-ray  and  electron  tubes,  the  development  of  quantum 
mechanics,  and  the  complete  theory  of  atomic  structure.  Today  It  can  be 
safely  claimed  that  gas  discharge  is  no  longer  a mystery:  no  unknown 

fundamental  mechanisms  are  believed  to  be  at  work  behind  the  phenomenon. 
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All  participating  physical  processes  are  accounted  for.  The  general  gas 
discharge  theory  can  be  formulated  in  terms  of  partial  differential  equations. 
Nevertheless,  it  Is  quite  another  matter  to  derive  particular  solutions  of 
these  equations  describing  particular  modes  of  the  discharge,  such  as  the 
corona.  As  the  equations  are  coupled  and  nonlinear,  one  has  here  the 
essential  trappings  of  a very  difficult  mathematical  problem. 

This  report  studies  various  solutions  of  the  basic  gas  discharge  equa- 
tions applicable  to  the  corona  on  a long  thin  wire.  Section  2 discusses  the 
observational  aspects  of  the  electrical  corona  and  the  physical  processes 
that  take  place  In  It.  A set  of  basic  equations  are  written  down.  Section 
3 presents  the  solutions  of  two  linear  boundary-value  problems  of  the  basic 
equations  to  determine  the  corona  breakdown  field  strengths  for  the  coaxial 
cylindrical  electrode  geometry,  both  for  quasi-static  and  for  microwave 
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excitations.  Section  4 solves  three  nonlinear  boundary-value  problems 
to  calculate  the  post-breakdown  d-c  corona  voltage-current  characteristic 
for  the  coaxial  electrode  configuration,  taking  successively  Into  account 
three  major  physical  effects.  Section  5 extends  the  d-c  corona  voltage- 
current  characteristic  calculations  to  the  line-over-plane  (or,  equivalently, 
parallel-line)  geometry.  The  results  are  applied  to  set  up  and  solve  non- 
linear transmission  line  equations  describing  the  d-c  corona  effect  on 
signal  propagation  along  a single  line  above  a conducting  ground  plane. 

Finally,  Section  6 discusses  the  limitations  of  the  present  results  and  the 
importance  of  some  as  yet  unsolved  problems. 

It  may  be  appropriate  to  state  at  the  outset  what  these  limitations  and 
unsolved  problems  are.  The  present  study  forms  the  preliminary  phase  of  an 
effort  to  determine  the  corona  effect  on  aircraft  long-wire  antennas.  Its 
goal  is  to  establish  the  basic  equations  governing  the  general  corona  phe- 
nomenon, and  to  understand  the  mathematical  nature  of  the  antenna  corona 
problem.  A set  of  boundary-value  problems  of  these  basic  equations  are  solved, 
and  the  solutions  compare  favorably  with  measurement.  However,  the  results 
are  valid  only  in  the  d-c  or  quasi-static  regime.  In  the  case  of  a long-wire 
antenna  excited  by  an  EMP  or  driven  by  an  LF  transmitter,  the  corona  problem 
is,  in  general,  a time-dependent  boundary-value  problem.  The  mathematical 
difficulty  involved  will  be  considerable.  Its  solution  must  be  relegated  to 
future  investigations.  More  discussions  on  this  point  will  be  found  below  in 
Section  6. 
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SECTION  2 

PHYSICS  OF  THE  CORONA 

2.1  Phenomenology 

The  electrical  corona  can  be  generated  by  the  following  apparatus 
arrangement  In  the  laboratory.  Take  a long  piece  of  thin  well-polished 
wire  and  encase  It  In  a large  concentric  cylindrical  metallic  shell. 

Apply  a static  voltage  difference  across  the  two  conductors,  with  the 
wire  positive.  As  the  voltage  Is  gradually  raised  and  approaches  a 
certain  critical  value,  the  electric  field  at  the  wire  becomes  very 
strong.  Substantial  Ionization  of  the  air  molecules  in  the  neighborhood 
takes  place.  When  the  critical  threshold  is  exceeded,  a luminous  bluish 
uniform  sheath  appears  on  the  wire  surface,  forming  the  corona  proper. 

Its  thickness  grows  with  Increasing  applied  voltage.  The  sheath  Is 
essentially  a highly  Ionized  plasma  with  positive  and  negative  Ions 
coexisting  in  comparable  concentrations.  The  luminosity  is  a tell-tale 
sign  of  Ion  recombination  and  neutralization  through  the  emission  of 
photons.  A shunt  current  can  be  detected  to  flow  from  the  wire  through 
the  corona  sheath  and  the  dark  space  beyond  Into  the  outer  cylinder. 

This  Indicates  a dielectric  breakdown  of  the  air. 

If  now  the  experiment  is  repeated  with  the  wire  negative,  the  corona 
Is  observed  to  assume  the  form  of  Isolated  reddish  tufts  cr  beads  more  or 
less  uniformly  spaced  along  the  wire.  The  breakdown  voltage  Is  also 
slightly  off  the  positive  corona  value.  The  sign  and  size  of  the  difference 
turn  out  to  depend  on  gas,  pressure  and  wire  material.  If  the  experiment 
Is  performed  with  an  alternating  voltage  difference,  the  a-c  corona  appears 
only  during  that  fraction  of  the  period  when  the  voltage  exceeds  the  critical 
value.  Its  visual  characteristics  exhibit  the  same  variation  with  wire 
polarity  as  for  the  d-c  corona.  Essentially  the  same  corona  observations 
can  be  made  with  the  parallel-wire  electrode  arrangement.  It  must  also  be 
remarked  that  If  Che  electrode  spacing  Is  very  small,  a spark  discharge  Is 
observed  and  the  corona  Is  not  formed. 

Because  of  the  possibility  of  substantial  corona  loss  in  high-voltage 
transmission,  Peek^  early  this  century  made  a series  of  careful  measure- 
ments on  long-wire  corona  characteristics.  From  the  data  he  derived  a set 


of  empirical  formulas  which  are  much  respected  even  today.  One  set  of 
formulas  gives  the  breakdown  electric  field  as  a function  of  the  wire  radius; 
another  set  gives  the  corona  loss  as  a function  of  the  applied  voltage. 

These  are  described  below. 

Breakdown  Field.  At  normal  temperature  and  pressure  (T  = 25“C, 

p = 760  mm.  Hg)  and  an  applied  frequency  of  60  Hz,  the  breakdown  or  critical 

electric  field  E on  the  wire  surface  of  radius  a at  which  the  corona 
c 

is  Initiated  is  given  by 


31  (l+  °^°^)kV/cm. , 


a In  cm. 


for  the  coaxial  geometry,  and  by 


- 29.8[H-  °'-^^)kV/cm. 
/& 


for  the  parallel-wire  geometry.  It  is  well  to  note  that  E^  depends  only 
on  the  wire  radius,  irrespective  of  the  location  of  the  other  electrode. 

For  large  a,  E^  approaches  a limit  of  about  30  kV/cm.  which  is  often  re- 
ferred to  as  the  dielectric  strength  of  air.  For  small  a,  E^  diverges 
rapidly.  The  numerical  constants  in  (1)  and  (2)  vary  with  atmospheric 
condition.  The  quantitative  dependence  has  been  investigated  by  Peek,  and  more 

3 

fully  by  Stephenson. 

Corona  Loss.  Peek  measured  the  power  loss  due  to  corona  and  obtained 
the  following  empirical  loss  formula  for  a single  line  above  ground: 

P « 12(f+25)^^  (V-V^)^  xio"^°  watts/meter  (3) 

where  f is  the  applied  frequency  in  hertz,  a and  h the  wire  radius  and  its 
height  above  ground,  V and  are  the  peak  values  of  the  applied  and  critical 

voltages  in  volts  relative  to  ground.  is  related  to  the  parallel-wire  break- 

down field  in  (2)  by 

\ . (A) 
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Due  to  the  difficulty  in  making  precise  high-voltage  loss  measurements, 
formula  (3)  is  admittedly  not  very  accurate,  particularly  for  V close  to 
V^.  The  quadratic  variation  of  P with  the  over-voltage  is  not  always 
corroborated  by  the  findings  of  other  investigators.  The  noteworthy 
feature  of  (3)  is  the  linear  dependence  of  the  loss  on  the  applied  fre- 
quency, indicating  a constant  loss  per  cycle. 

2.2  Fundamental  Processes  in  the  Corona 

The  one  basic  physical  process  which  underlies  the  corona  phenomenon 
is  the  collision  ionization  of  air  molecules  by  free  electrons.  The 
latter  are  always  present  in  air  to  some  extent  as  a result  of  cosmic 
ray  bombardment,  natural  radioactivity,  air  friction,  and  many  other 
causes.  In  an  intense  electric  field  such  as  existing  near  a highly 
charged  conductor,  these  mobile  electrons  are  accelerated  and  quickly  pick 
up  sufficient  energy  to  break  up  or  ionize  air  molecules  by  collision, 
producing  each  time  a positive  ion  and,  what  is  Important,  one  or  more 
daughter  electrons.  These  new  electrons  in  turn  can  initiate  collision 
ionizations  on  their  own.  When  conditions  are  right,  a chain  reaction 
takes  place.  A large  number  of  ions  and  electrons  are  produced  in  an 
Instant,  making  the  air  conducting  and  able  to  sustain  a discharge. 

The  criterion  which  permits  the  rapid  build-up  of  charged  particles 
in  a small  region  of  space  is  that  the  local  rate  of  particle  production 
by  ionization  exceed  the  rate  of  loss.  The  primary  particle  loss 
mechanisms  are  drift  and  diffusion.  In  an  external  electric  field, 
charged  particles  drift  along  the  field  lines  while  suffering  collisions 
with  one  another  and  with  neutral  particles.  Diffusion  arises  from  the 
microscopic  random  motion  of  these  colliding  particles,  and  is  related  to 
gas  pressure.  It  forces  particles  to  move  from  regions  of  high  particle 
concentration  to  those  of  lower. 

Besides  the  primary  production  and  loss  mechanisms  of  electron 
collisicn  ionization,  charged  particle  drift  and  diffusion,  some  secondary 
mechanisms  are  worth  mentioning.  Charged  particles  can  also  be  produced 
by  ion  collision  ionization,  photo-ionization,  and  electrode  emission. 
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They  are  also  lost  by  electrode  absorption,  recombination  (electrons  and  positive 
Ions  forming  neutral  molecules) , and  attachement  (electrons  attaching  themselves 
to  neutral  molecules  to  form  negative  Ions) . 

2.3  Mathematical  Formulation 

It  remains  to  translate  the  above  physical  Ideas  into  precise  mathematical 
formulas.  The  most  powerful  analytical  tool  available  to  date  for  achieving 
this  end  Is  the  Boltzmann  transport  equation,  which  describes  the  time  development 
of  a large  aggregate  of  particles  In  terms  of  their  fundamental  microscopic 
Interactions.  This  klnetlc-theoretlcal  approach  will  nevertheless  not  be 
adopted  In  the  following  calculations.  A more  coarse-grained  formulation  will 
be  employed.  It  will  be  a continuum  theory  based  on  a set  of  hydrodynamic 
equations,  which  are  the  first  few  moment  equations  of  the  Boltzmann  equation 
closed  by  truncation.  These  equations  contain  a number  of  transport  coefficients 
which  can  either  be  calculated  from  first  principles  In  kinetic  theory  or  simply 
taken  over  from  laboratory  measurements. 

The  chief  constituents  of  air  are  oxygen,  nitrogen  and  hydrogen.  The 
charged  particles  In  air  are  therefore  largely  made  up  of  their  positive  and 
negative  Ions  and  electrons.  In  the  following  an  Idealized  model  of  Ionized 
air  will  be  used  to  make  the  calculations  manageable.  Ionized  air  will  be 
taken  to  consist  of  one  representative  species  of  singly-charged  positive  Ions 
and  one  representative  species  of  singly-charged  negative  Ions.  They  move  under 
external  and  mutual  forces  over  a background  of  neutral  air  molecules.  This 
Idealization  Is  by  no  means  an  oversimplification,  since  experimentally  It  Is 
often  the  averaged  properties  of  Ions  of  one  sign  that  are  directly  measured. 
Essentially  one  deals  here  with  a two-charged-fluld  model  of  Ionized  air.  The 
neutral  molecules  actually  do  not  appear  directly  In  the  picture.  The 
field  variables  are  the  particle  number  d^^ltles  n^(£,t)  of  the  two  charged 
fluid  components,  and  the  total  electric  field  ^(£,t).  The  effect  of 
magnetic  fields  on  particle  motion  Is  negligible.  The  three  variables  are 
coupled  together  by  three  sets  of  field  equations  to  be  described  below. 

Equations  of  Motion.  These  equations  express  the  particle  current 
or  f?.ux  densities  X+(r,,t)  In  terms  of  the  number  densities  n^  and  the 


7 


electric  field  J|.  As  mentioned  previously,  the  motion  of  charged  par- 
ticles In  air  consists  of  a drift  and  a diffusion.  The  flux  density  has 
two  corresponding  components: 


~ —drift  ^—diffusion 
The  drift  flux  is  given  by 


(5) 


— drift 


± KnE 


(6) 


showing  a direct  proportionality  of  the  drift  velocity  to  the  field.  The 
physics  behind  this  formula  is  that.  In  a weakly-ionized  dense  gas  con- 
taining a large  amount  of  neutral  molecules,  a charged  particle  moving  in 
an  electric  field  does  not  experience  a sustained  acceleration.  As  soon 
as  enough  kinetic  energy  Is  accumulated,  it  Is  lost  through  Inelastic 
collisions  with  other  particles.  In  fact,  most  energy  Is  lost  In  ionizing 
collisions  with  neutral  molecules.  The  effect  of  these  collisions  on  the 
atomic  level  Is  the  appearance  of  a macroscopic  viscous  force  dragging  the 
particle's  movement.  The  drag  soon  counterbalances  the  electric  field, 
and  the  particle  henceforth  moves  on  with  a steady  velocity  proportional 
to  the  field.  The  constant  of  proportionality  is  the  ion  mobility  K in 
Eq.  (6).  The  mobility  is  always  defined  positive,  so  that  a ± sign  Is 
needed  in  (6)  to  account  for  the  oppositely  directed  motions  of  positive 
and  negative  charges  In  the  field.  Actually  Eq.  (6)  is  nothing  but  Ohm's 
law  In  disguise.  This  drag  phenomenon  has  an  analogue  In  parachute  jump. 
After  an  initial  gravitational  acceleration,  the  viscous  force  of  air  soon 
equalizes  gravity  and  the  free  fall  proceeds  at  a uniform  terminal  velocity. 

The  diffusive  flux  density  accounts  for  the  random  motion,  due  to 
random  atomic  collision,  that  Is  left  out  In  (6).  It  is  given  by 


f 


diffusion 


-7(Dn) 


(7) 


where  D is  the  diffusion  coefficient.  In  his  investigation  of  Brownian 
motion,  Einstein  showed  that,  at  thermal  equilibrium, 

k T 

_D_  _ V (8) 

K " e 
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where  e is  the  electronic  charge,  T the  equilibrium  temperature,  and 
kg  the  Boltzmann  constant.  The  effect  of  diffusion  is  to  even  out  par- 
ticle densities,  leveling  off  high  concentrations  and  filling  up  deple- 
tions. In  fact  the  phenomenon  is  mathematically  equivalent  to  heat  con- 
duction. 

In  summary,  the  flux  densities  for  positive  and  negative  ions  are 

(9) 

(10) 


i+  “ V+i"  * 

f _ - -K_n_E -7(D_n_). 


Poisson's  Equation.  The  total  electric  field  consists  of  an  external 
component  due  to  the  electrodes  and  an  internal  component  due  to  the  ions. 
It  Is  determined  by  the  Poisson  equation: 


7 -E  - (n^-nj  (11) 

o 

with  appropriate  boundary  conditions  at  the  electrodes.  The  permittivity 
of  air  is  taken  to  be  essentially  the  same  as  that  of  free  space. 

Equations  of  Continuity.  These  equations  express  particle  conserva- 
tion. They  balance  out  particle  gain  and  loss  mechanisms,  and  read 


3n 

-g^+7*  f^  « S(n^,n_,E), 

(12) 

3n_ 

Tt”'*’^*--  “ S(n^,n_,E). 

(13) 

where  S is  the  net  rate  of  particle  production  dependent  on  n^  and  E. 

The  right-hand  sides  of  (12)  and  (13)  are  the  same  since  positive  and  negative 
ions  are  produced  or  destroyed  at  the  same  rate  in  the  same  ionization  or 
recombination  processes.  The  functional  form  of  S for  electron  collision 
ionization  will  be  given  below. 

Equations  (11),  (12)  and  (13)  form  a set  of  coupled  partial  differential 
equations  in  n^  and  E,  with  defined  by  (9)  and  (10).  The  coefficients 

and  are  usually  regarded  constant,  with  values  taken  from  experiment. 
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These  equations  are  nonlinear.  The  nonllnearltles  occur  In  the  products 
n^E^  and  In  S.  Strictly  speaking,  the  equations  should  be  supplemented 
by  another  set  describing  energy  conservation.  As  It  turns  out,  corona 

4 

energy  dissipation  proceeds  mainly  by  way  of  ohmic  loss  which  Is  already 
Implicit  In  £q.  (6). 

Collision  Ionization.  When  an  electron  moves  through  a gas  under 
an  electric  field,  the  number  of  new  electrons  It  produces  by  collision 
Ionization  along  a unit  length  of  Its  trajectory  is  experimentally  found 
to  be  directly  proportional  to  the  local  electron  density; 

dn  “ an  dl  (14) 

As  is  well-known,  this  type  of  relation  produces  exponential  multiplication. 
The  constant  of  proportionality  a Is  known  as  Townsend's  first  Ionization 
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coefficient.  It  is  a nonlinear  function  of  the  local  electric  field,  the 
nature  of  the  gas,  and  the  gas  pressure.  For  a d-c  field  It  Is  given  to  a 
good  accuracy  by  the  empirical  formula: 

a « Ap  exp  (-Bp/E)  (15) 

where  p Is  the  gas  pressure  and  A and  B are  empirical  parameters.  For  air, 
Che  latter  are  found  to  be 

A > 15  Ion  palrs/cm.  x mm.  Hg  (16) 


B ■ 365  volts/cm.  x mm.  Hg 


Then  Che  Ion  production  rate  S due  to  collision  Ionization  Is  given  by 


S - a f 


In  the  following  sections  the  above  equations  will  be  solved  to  study 
various  aspects  of  the  corona  discharge  from  a long  wire.  It  Is  Co  be 
noted  that  essentially  the  same  set  of  equations  occurs  in  the  theory  of 
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transistor  p-n  junction  operation,  the  difference  being  in  the  relative 
magnitudes  of  the  coefficients  and  parameters  and  in  boundary  conditions. 
The  corona  breakdown  has  a transistor  analogue  known  as  the  avalanche 
breakdown  of  the  p-n  junction. 
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CALCULATION  OF  BREAKDOWN  FIELD 

This  section  contains  two  calculations  to  determine  the  corona  breakdown 
field  In  the  quasl-statlc  and  microwave  limits.  The  calculations  Involve  the 
solution  of  two  boundary-value  problems  for  the  Ion  equations  (12)  and  (13) . 

They  will  be  carried  out  for  the  coaxial  electrode  configuration  shown  In  Fig.  1. 

This  highly  symmetric  geometry  both  simplifies  the  mathematics  and  makes  possible 
comparison  with  measurements,  as  It  Is  a favorite  choice  of  the  experimentalists. 

This  type  of  calculation  Is  substantially  simpler  than  the  full-fledged 
corona  problem.  It  deals  mainly  with  the  pre-breakdown  regime  when  very  few 
Ions  are  formed.  The  electric  field  is  not  distorted  by  space  charge,  and 
nonllnearltles  do  not  occur.  Consequently  the  Poisson  equation  Is  not  needed. 

Furthermore,  variations  with  time  can  be  left  aside  If  one  Is  not  concerned  with 
the  temporal  development  of  the  breakdown.  The  problem  Is  mathematically  one 
of  determining  the  eigenvalues  of  linear  ordinary  differential  equations. 

3.1  Quasl-Statlc  Breakdown 

Suppose  the  wire  In  Fig.  1 Is  kept  at  a positive  static  potential  V 
relative  to  the  outer  cylinder.  Before  breakdown  there  are  very  few  Ions  in  the 
Intervening  air,  and  the  electric  field  is  given  by  the  familiar  expression 

E - E4„,  E(p)  - E(a)  i E(a)  - . (19) 

The  Ionization  coefficient  a In  Eq.  (15)  becomes  a function  of  space: 

o(p)  - Ap  exp[- . (20) 

Let  a free  electron  be  generated  by  some  means  at  the  surface  of  the  outer  electrode. 

It  will  travel  radially  Inwards  under  the  field,  and  on  the  way  produce  new 
electrons  by  Ionizing  collisions  with  air  molecules.  Upon  reaching  the  Inner 
wire  It  will  have,  according  to  (14) , multiplied  Into  an  avalanche  consisting 
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Figure  1.  Coaxial  electrode  configuration  for  corona  study. 


of  N electrons,  with 


N - exp[  o(p)dp] . 


At  the  same  time,  N-1  positive  Ions  are  produced  which  migrate  outwards  under 
the  field,  and  are  eventually  neutralized  and  absorbed  at  the  outer  cylinder. 
Individual  electron-initiated  avalanches  can  actually  be  observed  and  photographed 
In  the  laboratory.  The  avalanche  constitutes  a discharge  current. 

If  nothing  else  occurs,  a single  electron  can  only  trigger  a short  discharge 
pulse  this  way.  For  a continuous  discharge  to  take  place,  there  must  be  some 
mechanism  to  replenish  the  electron  supply  at  the  outer  cylinder,  which  In  this 
case  Is  the  cathode.  Townsend^  suggested  the  following  cathode  process:  upon 
Impact  on  the  cathode,  each  of  the  N-1  positive  Ions  has  on  the  average  a probability 
Y of  liberating  an  electron  from  the  cathode.  Therefore  a total  of  y(N-1) 
free  electrons  will  be  liberated.  It  Is  thus  clear  that,  starting  with  one  free 
electron  at  the  cathode,  one  ends  up,  after  Ionization  and  cathode  processes, 
with  y(N-1)  free  electrons  again  at  the  cathode.  If  y(N-1)  ■ 1,  one  ends  up 
precisely  with  what  one  begins.  The  cycle  can  repeat  Itself  endlessly,  and  the 
discharge  becomes  self-sustaining.  A continuous  discharge  current  flows 
between  the  electrodes.  This  Is  the  point  when  the  corona  starts  to  appear. 

By  (21) , the  condition  for  this  occurrence  Is 


fb 

(exp  ad  p - l]  ■ 1 . 

J A 


It  Is  known  as  Townsend's  breakdown  criterion.  The  coefficient  y Is  dependent 
on  gas,  cathode  material,  and  other  parameters. 

The  condition  (22),  derived  above  by  simple  physical  arguments,  can  also 
be  obtained  formally  by  solving  a boundary-value  problem  of  the  Ion  equations  (12) 
and  (13).  The  time-independent  form  of  these  equations  for  the  coaxial  geometry 


-«(p)|f.(p)l 
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^ 4- [»'-'»>]  - «(p)lf.(p)l- 


(24) 


It  Is  convenient  to  Introduce  the  total  positive  and  negative  Ion  fluxes  F^(p) 
per  unit  length  of  wire  such  that 

F.(o)  F (p) 

■ - iiT  ■ ^ “■  <• 


Then  (23)  and  (24)  become 


^ ^^.(p)  “ oF_(p), 


- ^ F_(p)  - aF_(p), 


The  solution  Is 


F (p)  - F (b)  exp 


F_j^(p)  + F_(p)  » constant.  (28) 

The  constants  of  Integration  are  determined  by  boundary  conditions  at  the  electrodes. 
At  the  cathode  p ■ b,  one  has  YF_j_(b)  • F_(b)  as  discussed  above.  At  the  anode 
p * a,  namely  the  wire,  the  positive  Ion  flux  must  be  0 since  these  Ions  are 
‘generated  In  air  and  migrate  towards  the  cathode.  Hence  one  has  F^(a)  > 0. 
Substitution  of  (27)  and  (28)  Into  the  boundary  conditions  yields  the  following 
pair  of  homogeneous  linear  algebraic  equations: 


F_(a)  ■ F_(b)  exp 


(l+  ^)F_(b)  - F_(a). 


These  equations  have  the  trivial  solution  (namely  zero  Ion  current)  unless  the 
determinant  vanishes.  This  condition  Is  easily  seen  to  be  Identical  to  (22). 
In  other  words,  the  solution  describing  the  breakdown  of  the  air  gap  between 
the  electrodes  Is  the  elgen-solutlon  of  Eqs.  (23)  and  (24). 
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At  breakdown  the  field  at  the  wire  E(a)  is  by  definition  the  breakdown 
field  E^.  Substitution  of  (20)  into  (22)  yields 


This  relates  E^  to  the  geometric,  atmospheric,  and  other  parameters  of  the 
problem.  The  solution  showing  the  variation  of  E^  with  a is  plotted  in  Fig.  2 
for  A and  B given  by  (16)  and  (17),  p “ 760  mm.Hg,  y * 0.025  for  a copper 
cathode  in  air^,  and  b/a  > 100  as  a typical  radii  ratio  for  the  electrodes. 
Actually  the  result  is  totally  insensitive  to  this  last  ratio,  as  the  second 
exponential  in  (31)  is  entirely  negligible  for  b only  a few  times  larger 
than  a.  The  breakdown  data  of  Peek  taken  for  concentric  copper  cylinders  at 
60  Hz'  are  also  plotted  for  comparison.  The  data  points  follow  the  theoretical 
curve  closely,  making  the  same  dramatic  upturn  for  a <0.2  cm.  The  measured 
values  are  about  10%  lower  than  predicted.  In  view  of  the  crudeness  of  the  theory, 
which  accounts  only  for  ionization  and  cathode  processes,  the  agreement  between 
experiment  and  theory  can  be  considered  satisfactory.  Much  better  agreement 
can  certainly  be  achieved  by  adjusting  the  parameters  in  (31)  to  fit  the  data. 

The  above  calculations  are  performed  for  a positive  wire.  It  can  be  shown 
that  exactly  the  same  breakdown  criterion  is  obtained  for  a negative  wire.  This 
theory  therefore  cannot  account  for  the  small  dependence  of  the  breakdown  field 
on  wire  polarity. 


3.2  Microwave  Breakdown 

In  an  applied  electric  field  at  microwave  frequency,  an  ion  moves  quite 
differently  than  in  a d-c  field.  The  field  reverses  its  direction  every  half -cycle, 
and  the  ion's  motion  is  predominantly  oscillatory.  There  is  very  little  drifting 
motion,  the  mobility  being  inversely  proportional  to  the  applied  frequency. 

In  this  limit  the  equation  of  continuity  (13)  for  electrons  becomes 

3n  » 

-7^(D_n_)  - S.  (32) 
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Lgure  2.  Breakdown  electric  field  in  air  at  normal  temperature  and  pressure  as  a function 

of  wire  radius  a.  The  crosses  are  data  of  Peek  (ref.  1)  taken  at  60  Hz. 


S.C.  Brown  and  coworkers  have  applied  this  equation  to  calculate  the  breakdown 
field  in  microwave  cavities,  and  verified  the  results  by  experiment.  In  their  work 
they  took  for  S,  the  rate  of  electron  production  by  collision  ionization,  the 
following  form  equivalent  to  (18) : 


S « v^(E,a))n^, 


where  v^,  is  the  ionization  frequency,  in  general  field-  and  frequency-dependent. 


With  the  introduction  of  a function  ^ such  that 


= D n , 


Brown  obtained  from  (32)  the  time-dependent  diffusion  equation 


+ v^ip. 


Before  breakdown,  there  are  very  few  electrons  and  ions  in  air.  The  electric 
field  E is  essentially  the  same  as  the  externally  applied  cavity  field.  Its 
spatial  dependence  is  experimentally  controllable.  Then  Eq.  (35)  is  linear  and 
homogeneous,  and  can  be  studied  with  the  method  of  separation  of  variables  by 
setting 


t|j(r,t)  - I tiO. 


Substitution  of  (36)  into  (35)  yields  the  equation 


(D  V + V.  + \ )^  =.  0 

— 1 m m 


for  the  eigenfunctions  iji^.  These  and  the  corresponding  eigenvalues  are 

determined  by  the  geometry  of  the  microwave  cavity  and  the  boundary  conditions. 

Brown  assumed  that  the  cavity  walls  were  perfectly  absorbent  to  impinging  electrons. 
In  a diffusion  process  this  means  that  the  electron  density  will  stay  very  low  in 
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the  vicinity  of  the  wall  surface,  and  will  actually  extrapolate  to  zero  a short 
distance  behind  the  wall.  Therefore  the  boundary  condition  was  taken  to  be 


0. 


at  cavity  walls. 


(38) 


This  uniquely  determines  the  eigenfunctions  and  the  eigenvalues. 

For  a general  microwave  cavity,  the  breakdown  condition  can  be  expressed  in 
terms  of  the  stability  of  the  solution  (36).  Before  breakdown,  0 is  stable. 

This  requires  that  all  the  eigenvalues  be  positive.  It  is  clear  that  as 

the  breakdown  threshold  is  approached,  the  eigenvalue  spectrum  is  lowered;  and 
the  lowest  eigenvalue  (say,  X^)  approaches  0.  For  X^  only  slightly  negative, 
in  (36)  is  unstable  and  grows  exponentially  with  time.  The  breakdown  criterion 
is  therefore  that  the  lowest  eigenvalue  be  0 or,  equivalently,  that  the  equation 


(D_V^  + = 0 (39) 

have  a non-trlvlal  solution. 
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Berlin  and  Brown  considered  Eq,  (39)  for  a coaxial  microwave  cavity  of 
finite  length.  The  cavity  was  excited  by  a definite  cavity  eigenmode,  with  definite 
spatial  dependence  of  the  electric  field.  The  breakdown  field  intensity  was 
measured  for  the  elgenfrequency.  Separately,  experiments  were  done  to  determine 
the  field-  and  frequency-dependence  of  the  coefficients  D_  and  v^.  The 
data  were  fed  into  (39)  in  the  form  of  empirical  formulas.  The  imposition  of  the 
non- triviality  condition  on  the  solution  led  to  a transcendental  equation  for  the 
breakdown  field.  Ref.  9 should  be  consulted  for  calculatlonal  details.  The 
calculated  breakdown  field  was  found  to  be  in  excellent  agreement  with  experiment, 
making  this  theory  perhaps  the  most  successful  breakdown  calculation  to  date. 

The  above  calculation  of  the  microwave  breakdown  field  is  analogous  to  the 
theoretical  determination  of  the  critical  mass  of  a bulk  of  uranium.  Neutrons 
in  a uranium  sample  diffuse  freely,  and  are  scattered,  captured  and  emitted  by 
the  uranium  nuclei.  They  Initiate  nuclear  reactions  with  the  nuclei,  resulting 
in  the  release  of  huge  amounts  of  energy  and  daughter  neutrons,  in  much  the  same 
way  as  collision  ionization  in  a gas  releases  new  electrons.  The  neutron  density 
satisfies  a diffusion  equation  similar  to  (35).  For  a uranium  sphere  with  boundary 
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condition  (38),  the  lowest  eigenvalue  turns  negative  when  the  sphere  radius 
exceeds  a certain  critical  value.  The  solution  becomes  unstable,  and  the 
neutrons  multiply  exponentially.  A nuclear  chain  reaction  is  set  off  with  explo' 
sive  consequences. 
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CALCULATION  OF  CORONA  VOLTAGE-CURRENT  CHARACTERISTIC 

Consider  the  coaxial  electrode  configuration  depicted  in  Fig.  1. 

When  the  voltage  difference  applied  across  the  gap  exceeds  the  critical 
value,  the  corona  appears  on  the  wire  surface.  The  insulation  provided 
by  the  air  ceases,  and  a current  flows  between  the  wire  and  the  outer 
cylinder.  For  practical  applications  it  is  important  to  know  the  depen- 
dence of  this  corona  current  on  the  applied  voltage,  or  the  voltage- 
current  characteristic.  This  section  is  devoted  to  its  calculation  in 
the  quasi-static  limit. 

As  mentioned  previously,  the  corona  sheath,  which  clings  to  the  wire, 
is  a plasma  in  which  electrons  and  positive  ions  coexist  in  comparable 
concentrations.  For  a positive  corona,  the  electrons  are  accelerated 
towards  the  wire  and  the  positive  ions  are  repelled  towards  the  cylinder. 
Thus  one  expects  to  find  a layer  of  electrons  bound  to  the  wire  and  a 
stream  of  positive  ions  flowing  from  the  wire  to  the  cylinder.  The  corona 
sheath  indicates  the  extent  of  the  electron  layer,  and  the  corona  current 
is  mainly  carried  by  positive  ions.  For  a negative  corona,  one  would 
expect  the  electrons  and  the  positive  ions  to  exchange  roles.  There  is, 
however,  a slight  complication:  as  soon  as  the  electrons  leave  the  sheath 
and  move  towards  the  cylinder,  they  attach  themselves  to  neutral  molecules 
to  form  negative  ions.  Therefore,  for  the  negative  corona,  the  corona 
current  is  carried  by  electrons  inside  the  corona  sheath  and  by  negative 
ions  outside. 

An  Important  conclusion  can  be  drawn  from  the  above  considerations: 
in  the  quasi-static  limit  the  corona  current  is  predominantly  carried  by 
charges  of  the  same  sign  as  the  wire.  As  an  approximation,  one  neglects 
the  layer  of  bound  charges  in  the  corona  current  calculations. 

4.1  Effect  of  Ion  Mobility 

Just  as  in  metallic  conduction,  the  chief  determining  factor  in  gaseous 
conduction  is  the  mobility  of  the  carriers.  Taking  only  this  effect  into 
account,  one  can  calculate  the  d-c  corona  current  for,  say,  the  positive 
corona  from  particular  cases  of  Eqs.  (9),  (11)  and  (12): 
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(40) 
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(41) 

^ * 1+*  0 

(42) 

These  equations  describe  a time- Independent  conduction  model  In  which 
positive  Ions  Interact  self-conslstently  with  the  electric  field,  but 
are  neither  produced  nor  destroyed  by  Ionization  and  recombination  pro- 
cesses. This  model  appears  to  have  been  first  considered  by  Townsend^^ 
For  a cyllndrlcally  symmetric  geometry,  the  thre®  equations  become 
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The  Integration  of  (45)  yields 


f 


+ 


1 

2ireo 


(46) 


where  1 Is  a constant  of  Integration  and  has  the  meaning  of  the  corona 
current  per  unit  length  of  wire.  Elimination  of  f^  between  (43)  and  (46) 
gives 


en. 


2itK  pE 


(47) 


Substitution  of  this  Into  the  right-hand  side  of  (44)  results  In  the 
nonlinear  differential  equation: 


ire  K . 
o + 


(48) 
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Integration  from  a to  p leads  to 


2 2 2 2 1 2 2 
p^E  (p)  - aV(a)  + ^ (p^-a^).  (49) 

Solving  for  E(p)  and  Integrating  across  the  gap  from  a to  b,  one  obtains 
the  voltage  difference 


V 


aE(a)  V (X^-1>^+1  - 1 - In 


ML 


X(./I^  + l) 


(50) 


where 


X 


*0  ■ 


(51) 


Eq.  (50)  la  a relation  between  the  applied  voltage  V,  the  corona  current  per  unit 
length  of  wire  1,  and  the  electric  field  on  the  wire  E(a).  In  order  to  arrive 
at  a unique  voltage-current  characteristic,  the  last  quantity  must  be  determined 
Independently.  For  this  purpose  one  Imposes  the  corona  boundary  condition: 

E(a)  ■ E^,  after  breakdown.  (52) 


That  Is  to  say,  the  field  on  the  wire  does  not  change  x^th  Increased  applied  voltage, 
but  persists  at  the  breakdown  value  E^.  The  system  adjusts  Itself  to  the  Increased 
voltage  by  emitting  an  appropriate  corona  current.  With  this  condition,  Eq.  (50) 
becomes 


r rr 


l)(f+l  - 1 + (/i^ - l)je,nx  - to 


(53) 


where  Is  the  breakdown  voltage  given  by 


V - aE  lln(b/a) 
c c 


(5A) 


for  the  coaxial  configuration. 
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The  corona  voltage-current  characteristic  (53)  is  plotted  in  Fig.  3 

2 

for  the  case  100.  Typically  E^'\<  30  kV  /cm.,  'v  1 cm.  /volt  sec. 

From  this  one  has  i^ 0.1  A/m.  Fig.  3 shows  that  a typical  corona  current 
is  of  the  order  of  1 mA/m.  The  curve  i versus  V has  a parabolic  shape, 
and  can  be  fitted  closely  by  the  simple  form 


i - icV  (V-V  ) 
c 


(55) 


This  was  first  obtained  by  Townsend,  and  has  been  periodically  verified 

by  experimentalists.  Fig.  4 shows  some  recent  data  reported  by  Albrecht, 

Wagner  and  Bloss  for  the  positive  corona.  The  measurements  were  done  in 

pure  nitrogen  for  the  coaxial  electrode  geometry  at  1 atmosphere.  The 

current  carriers  are  ions.  The  corona  voltage-current  characteristic 

is  plotted  with  the  current  on  a logarithmic  scale,  and  compared  with  the 

theoretical  expression  (53).  Excellent  agreement  is  indicated. 

Fig.  5 shows  the  radial  variations  of  the  electric  field  E(p),  the 

potential  V(p),  and  the  ion  density  n^(p)  for  the  case  V ■ 2V^.  From 

Fig.  3,  the  corresponding  corona  current  is  equal  to  0.005  1 . The  above 

o 

three  quantities  are  normalized  to  their  respective  boundary  values  E^,  V 
and  n^,  respectively,  on  the  wire  where 


n 

o 


e E 
o c 

ea 


(56) 


The  negative  corona  voltage-current  characteristic  can  be  obtained  in 
a similar  fashion.  The  aforementioned  difference  in  the  nature  of  the 
negative  charge  carriers  inside  the  corona  sheath  from  that  outside  can 
be  accounted  for  by  an  effective  overall  negative  ion  mobility. 


4.2  Effect  of  Ion  Diffusion 

In  Eq.  (40),  the  part  of  the  particle  flux  due  to  diffusion  is  neg- 
lected. This  subsection  studies  the  voltage-current  characteristic  when 
the  diffusion  flux  is  retained.  The  total  flux  is  therefore  given  by  Eq. 
(9).  In  place  of  Eqs.  (43),  (44)  and  (45),  one  now- has 
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Nomalized  electric  field  E(p),  Ion  number  density  n(p),  and  potential  V(p)  across 
coaxial  electrodes  in  corona  with  applied  voltage  difference  V - 2V  . 


f+-  W-D+d^V 


where  for  simplicity  the  diffusion  coefficient  D_j_  is  considered  constants 
The  solution  of  (59)  is  again  (46).  Substitution  of  (46)  into  (57)  gives 


-D^  ^ n^  + K . n . E - -r^  . 
+ dp  + + + 21160 


Elimination  of  n_j^  between  (58)  and  (60)  leads  to 


where 


d u . du  . , du  2k 

^ , 2 dp  dp  A 

dp 


u « pE,  A 


2 

2D  ' ^ " 2 ’ 

8ire  DT 
o + 


Eq.  (61)  has  been  studied  by  Borgnis'*’^  in  his  investigation  on  electrolyte 
conduction.  The  following  analysis  is  largely  patterned  after  his.  The 
equation  is  equivalent  to 


d r du  . , . .,2,  2k 
_ [-P  _ + 2u+Au  1 - — P 


and  hence  a first  integral  is  the  Rlccatl  equation 


-p  ^ + 2u  + Au^  ■ ^ p ^ + C . 
dp  A 
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With  the  nonlinear  transformation 


A dp 


(64)  goes  over  to  a linear  equation  for  the  new  unknown  v: 


2 d V dv  /■,  2 2 . - 

P — 2 “ P ^ P +AC)v  ■ 0. 
dp 


This  Is  a modified  Bessel  equation,  and  the  solution  Is 


V - p[C'I^(kp)  +C"K^(kp)], 


where 


The  electric  field  Is  obtainable  from  (62)  and  (65)  as 


E ■ - Y 

A dp 


The  solution  contains  three  unknown  constants  C,  C"  and  v.  Since 
only  Inv  appears  In  (69),  It  Is  clear  that  only  the  ratio  C'/C"  matters. 
This  still  leaves  one  with  two  constants — one  more  than  E(a)  In  the 
previous  subsection.  The  reason  Is  that  the  Introduction  of  the  diffusion 
term  In  Eq.  (60)  Increases  the  order  of  the  differential  equation  by  1. 

An  extra  boundary  condition  pertinent  to  the  diffusion  process  must  be 
stipulated.  Instead  of  seeking  such  an  extra  diffusion  condition,  which 
Is  often  difficult  to  establish,  a different  approach  will  be  adopted  here. 

It  will  be  assumed  that  the  correction  due  to  diffusion  Is  small  (as  Is  amply 
Indicated  by  the  good  agreement  of  the  results  In  Subsection  4.1  with 


experiment),  so  that.  In  the  limit  0,  the  solution  should  tend  con- 

tinuously to  the  previous  diffusion-free  expression.  It  can  be  shown  that 
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for  a positive  corona,  Che  function  I ^ In  (67)  should  be  eliminated.  This 
essentially  gets  rid  of  the  extra  constant  of  Integration.  Therefore  one 
obtains  for  the  electric  field 


E - - ln[pK^(kp)].  (7( 

Integration  over  p from  a to  b yields  the  voltage  difference 


The  surviving  unknown  constant  v is  determined  from  the  corona  boundary 
condition  (52)  which  by  (70)  reads 


kaK' (ka) 

V 

K (ka) 

V ' 


AaE  +1. 
c 


(77 


Eqs.  (71)  and  (72)  together  constitute  the  corona  voltage-current  characteristic 
with  diffusion. 


It  can  be  shown  that.  In  the  limit  -►  0,  all  three  constants  A,  k 


and  V diverge.  One  can  therefore  use  Che  Debye  as3rmptotlc  formulas  for 
large  order 


-VC 


1/4 


-vt 


where 


C ■vl+2  + In- 


1+ Vl+2^ 


Substitution  of  (73)  Into  (72)  yields  the  simple  solution 

V ■ 


[AaE^+l]^  -k^a^ 


(7: 


(7- 


(7 


In  terms  of  the  reduced  corona  current  introduced  in  (51) , one  finds 


(ka)^  - (AaE  )^(ji  . 
c 


(7 
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Hence 


V - AaE  V(l+5)^-(j>,  5 - • 

c 


(77) 


In  the  Debye  asymptotic  limit,  the  voltage  difference  (71)  reduces  to 


V - j rVv2+(itb)^  - '^v^+(ka)^  - ^ ^ ^ 

^ vfvv^+(ka)^  - 


(78) 


From  Eqs.  (75),  (76)  and  (77),  It  can  be  seen  that  (78)  differs  from  (50)  only 
by  the  presence  of  the  quantity  6 in  (77).  By  (62)  and  (77),  one  obtains 


6 


2D, 


K^aE 
-♦*  c 


(79) 


2 2 - 

In  air  D^''- 0.028  cm.  /sec.,  K^''^1.36  cm.  /volt  sec.,  and  so  A * K^/2D^ 24/volt. 

It  then  turns  out  that  6 ''-10  which  is  entirely  negligible  compared  to  1 

in  (77).  Therefore  it  is  concluded  that  the  effect  of  diffusion  is  infinitesimal 

in  the  present  case  compared  to  that  of  mobility,  although  the  former  actually 

overtakes  the  latter  in  the  microwave  regime. 

The  effect  of  diffusion  on  the  voltage-current  characteristic  has  also 

13 

been  studied  by  Chekmarev,  but  with  a boundary  condition  other  than  (52) . 


4.3  Effect  of  Ionization 

The  two  calculations  of  the  corona  voltage-current  characteristic  presented 
above  are  based  on  the  assumption  that  the  current  is  carried  by  charges  of 
one  sign  emitted  by  the  wire.  This  is  called  the  unipolar  theory.  In  reality, 
since  the  corona  sheath  contains  charges  of  both  signs  created  by  ionization, 
both  types  of  charges  must  participate  in  conduction.  For  example,  in  a 
positive  corona,  some  of  the  electrons  in  the  sheath  will  be  attracted  to  the 
positive  wire  and  migrate  towards  it,  forming  a small  electron  current.  A 
calculation  which  takes  this  aspect  into  account  is  a bipolar  theory. 
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There  Is  as  yet  ao  reliable  quantitative  theory  for  the  state  of  affairs 
Inside  the  corona  sheath.  The  precise  Ionization  and  recombination  rates  and 
the  detailed  field  and  charge  distributions  Inside  the  sheath  are  difficult 
to  evaluate.  In  the  following  calculations  the  effect  of  Ionization  will  be 
represented  by  a constant  overall  Ionization  rate.  Therefore,  although  the 
boundary-value  problem  so  formulated  can  be  solved  exactly,  the  results  can  be 
expected  to  have  only  qualitative  and  methodological  values.  With  the  total 
neglect  of  diffusion  justifiable  from  the  results  of  the  last  subsection,  the 
basic  equations  are 


1+  “ 

1.  * 

• 

(80) 

7 . E - ~ 
— e 

(n^-  n_). 

(81) 

o 

7-1+  * s. 

7 • f _ - S 

(82) 

\ where  S Is  a constant  Ionization  rate 

. In  the  case  of 

axial  symmetry. 

these 

1 equations  become 

1 

f_  « -K_n_E 

(83) 

) 

f-  (n^-  n_) 

(84) 

: 

o 

! 

ii'-v  -5- 

(85) 

The  solutions  of 

(83)  and  (85)  are 

„ ..  1 2 

i 

(86) 

pK_^n_^E  “ 2 Sp  + C, 

K^n  E + K n E - 
+ + - - 

2irep 

where  C and  1 are  constants  of  Integration,  and  1 again  has  the  meaning  of 
the  corona  current  per  unit  length  of  wire.  For  a positive  corona,  the  negative 
charges  are  created  by  ionization  In  the  space  between  the  electrodes  and  migrate 
towards  the  wire.  The  density  of  negative  charges  n_  at  the  outer  cylinder 
p ■ b can  be  set  equal  to  0.  This  boundary  condition  determines  C.  Therefore 
the  solution  (86)  becomes 
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pK^n^E  - ^ - I (b^-P^).  pK_n_E  - | (b^-p^) . (87) 

Solving  for  n^  and  n and  substituting  into  (84),  one  obtains  the  differential 
equation 


_ d , i eS  f 1 1 •v  ,.2  2. 

dp  " 2ire  K,  ^ “ 2e  ^K.  "^K  ^ 


(88) 


'o  + 


o + 


Integrating  from  a to  p,  one  obtains 


(PE)‘  - (aE^)-+  K 


J O')  ®S(K  + K)  999  1 iA 


(P  -a")  - 


o + 


2tKK  rb^(P  - I(P  -)]  (89) 

o + - 


where  one  has  imposed  the  corona  boundary  condition  (52),  namely,  that  the 
electric  field  on  the  wire  surface  is  maintained  at  the  critical  value  E 

c 

after  breakdown.  Solving  (89)  for  E,  and  integrating  over  p from  a to  b, 
one  arrives  at  the  following  voltage-current  characteristic 


rb 


V =»  aE 


^ Vl  + 


- l]  - SW(p) 


where 


i/i  as  in  (51),  and 
o 


W(p) 


'«+*  >'->  ,.2,  2 2.  1 , 4 4. 


2 2 

20^a^E  K.K 
o c + — 


[b‘(p"-a‘)  - j (P  -a")] 


(90) 


(91) 


For  S ■ 0,  Eq.  (90)  reduces  to  the  lonizationless  limit  (53).  For  nonzero  S, 
the  integral  can  also  be  worked  out  analytically.  But  the  result  will  be  very 
complicated.  Since  ionization  effects  are  only  very  approximately  incorporated 
in  (85),  an  exhaustive  calculation  will  be  quite  out  of  place  here.  There  is, 
at  any  rate,  no  dependable  method  to  assign  an  effective  value  to  S.  A 
qualitative  examination  of  (90)  suffices.  It  is  observed  that  the  function  W(p) 
is  positive  inside  the  range  of  integration.  The  effect  of  ionization  is  therefore 
to  decrease  the  expression  under  the  square  root.  This  decrease  must  be  compensated 
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by  an  Increase  in  (j)  (or  the  corona  current  i)  if  the  voltage  difference 
V is  to  remain  the  same.  One  therefore  arrives  at  the  intuitively  obvious 
conclusion  that  ionization  augments  the  corona  current. 
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SECTION  5 

APPLICATION  TO  A SINGLE  LINE  ABOVE  GROUND 

In  the  previous  sections,  a number  of  boundary-value  problems  in  the 
study  of  the  electrical  corona  on  a long  wire  were  solved  for  the  coaxial 
electrode  configuration  illustrated  in  Fig.  1.  This  geometry  was  favored 
because  its  high  degree  of  symmetry  facilitates  the  mathematics,  and  because 
experimental  data  for  it  are  available  for  comparison  with  calculation.  In 
practice  there  exists  another  configuration  of  great  importance,  consisting  of 
a long  wire  parallel  to  the  surface  of  a ground  plane.  This  is,  for  example, 
the  geometry  of  the  high-voltage  power  line.  As  is  well  known,  if  the  ground 
plane  is  perfectly  conducting,  this  geometry  is  also  equivalent  to  that  of  two 
parallel  wires  operating  in  the  push-pull  mode  and  far  away  from  other  conductors. 

This  section  calculates  the  corona  voltage-current  characteristic  for 
the  wire-and-plane  configuration.  The  results  are  applied  to  study  the  corona 
effect  on  large  signal  propagation  along  a semi-infinite  wire  above  ground. 

5.1  Corona  Voltage-Current  Characteristic 

Consider  a long  straight  \d.re  of  radius  a whose  axis  lies  at  a height 
h above  a perfectly  conducting  ground  plane,  as  depicted  in  Fig.  6.  Suppose 
the  wire  is  charged  to  a potential  V relative  to  ground.  When  V exceeds  the 
corona  onset  value  V^,  a current  will  flow  through  the  air  between  the  wire 
and  the  ground.  The  critical  potential  V^  is  related  to  the  breakdown  field 
E^  through  Eq.  (4).  As  has  been  discussed  in  Section  3,  E^  is  a function  of 
a and  not  of  h.  It  is  for  all  practical  purposes  the  same  as  for  the  coaxial 
geometry. 

The  coordinate  system  most  suitable  for  studying  the  wlre-and-ground 
configuration  in  Fig.  6 is  that  of  the  bipolar  coordinates  (C,ri,z).  The  relations 
between  (x,y)  and  (CtH)  are 

(x  - d cot  n)  ^ + y^  ■ d^  csc^  n,  (92) 

x^  + (y  - d coth  5)^  - d^  csch^  5-  (93) 

The  parameter  d is  a positive  constant.  The  range  of  C is  from  to  ®, 


and  that  for  n from  0 to  2w.  Then  these  relations  describe  two  families 

of  orthogonally  Intersecting  circles  on  the  x-y  plane  as  shown  In  Fig.  7. 

Eq.  (92)  represents  a family  of  circles  with  centers  on  the  x-axls  and  passing 

through  two  common  points  y > td  known  as  the  foci.  Each  circle  Is  labeled  by 

a value  of  n.  Eq.  (93)  represents  the  orthogonal  net.  Each  of  Its  circles  Is 

labeled  by  a value  of  C.  The  upper  focus  Is  the  limit  5 and  the  x-axls 

the  limit  ^-^O.  The  circles  with  C>  0 are  below  the  x-axls.  Note  that  n ” ir 

for  X “ 0,  “d  < y d,  that  Is,  along  the  Interfocal  line. 

The  wire  In  Fig.  6 can  be  represented  by  one  of  the  circles  (93).  Let 

It  be  labeled  by  a value  C < 0.  It  can  be  shown  that 

o 


E --coah'^0. 


Ju2  2 
»h  - a . 


The  coordinate  n measures  the  angle  around  the  wire. 

In  the  following  corona  current  calculations,  only  the  effect  of  Ion 
mobility  will  be  considered.  The  basic  equations  are  the  same  as  Eqs.  (40), 
(41)  and  (42)  In  Subsection  4.1: 


7 . E - — n^, 
“ ^o 

^ " 0* 


(97) 


If  the  wire  Is  thin,  the  state  of  affairs  close  to  the  wire  may  be  considered 
axially  symmetric,  so  that  the  family  of  C-labeled  circles  are  equlpotentlal 
surfaces.  The  corona  current  flows  from  wire  to  ground  along  circles  of  constant 
n In  Fig.  7.  This  means  that  and  ^ have  only  5-components.  Writing 

out  the  above  equations  In  bipolar  coordinates,  one  obtains 


2 

(cosh  5 - cos  n)  3 f d -■>  _ e 

,2  35  ^cosh  5 - cos  n ’ 
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(100) 


r*^ 


2 

(cosh  C - cos  n)  3 f d f 1 - 0 

,2  35  ^cosh  5 ~ cos  n 


The  solution  of  (100)  Is 


C(n) 


(101) 


where  C(n)  depends  only  on  n*  For  a thln-wlre  the  corona  current  flowing  out 
from  the  wire  surface  is  isotropic.  This  means  that 


" 2Trea’  ^ " ^o 


(102) 


where,  as  before,  1 is  the  corona  current  per  unit  length  of  wire.  Under  this 
boundary  condition,  Eq.  (101)  becomes 


f m 1 cosh  5 - cos  n 
+ 2irea  cosh  5^  - cos  n’ 


Solving  for  n^  from  (98)  and  (103)  and  substituting  in  (99) , one  obtains 


- G 

ar  V 


35  '■cosh  5 ■ cos  n 


e) 


ire  a(cosh5  -cosn)  'ccsh  5 - cos  n 
o + o 


2 

) . 


Integration  from  5^  to  5 yields 


'cosh  5 ~ cos  n 


e)  - (; 


cosh  5 ^ cos  n 
o 


E ] 
.3 


JL 


ireK,  a(cosh5  -cosn)  J,  / u i-  \2’ 

o + o ^ -cosn) 


where  one  has  already  Imposed  the  corona  boundary  condition  (52),  that  is. 


E - E at  5 - C . 
c o 


(103) 


(104) 


(105) 
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In  bipolar  coordinates,  the  ^-component  of  the  field  Is  related  to  the 
potential  through 


E - - 


cosh  g - cos  n 3V 


3C’ 


(106) 


and  so  the  left-hand  side  of  (105)  Is  simply  the  square  of  3V/3?;.  The  voltage 
difference  can  be  obtained  by  Integrating  this  last  derivative  from  the  wire 
to  the  groimd.  The  path  of  Integration  Is  most  conveniently  chosen  as  the 
Interfocal  line  n “ v.  With  this  choice  the  Integral  In  (105)  becomes 


dC 


5^  (coshe  +1)‘ 


dC 


C 4cosh^(C/2) 
o 


j ^tanh(|-)  - y tanh' 


(107) 


Also,  by  (94), 


cosh  5+1 
o 


JEi 

IM-a  • 


(108) 


With  these  simplifications,  the  voltage  difference  between  wire  and  ground 
becomes 


^ r s 


(109) 


where  ■ ^^^o 


X(5)  - tanh(|-)  - j tanh^(-|-) 


(110) 


Just  at  breakdown,  V - V and  4>  ■ 0.  Eq.  (109)  gives 

c 


V - -C  aE  , 

c o c “h+a 


(111) 
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(112) 


V 

c 


aE 

c 


cosh"^(j) 


which  is  essentially  (4)  for  h » a. 

The  wlre-to-ground  corona  voltage-current  characteristic  (109)  is  computed 

for  the  case  h/a  ■ 100  and  plotted  in  Fig.  8.  It  has  the  same  parabolic 

appearance  as  that  in  Fig.  3 for  the  coaxial  geometry,  and  can  also  be  fitted 

by  the  simple  expression  (55).  In  fact,  it  is  obvious  from  (109)  that  for  a 

2 

larger  current  (it)  ■♦■“),  one  has  i^  V . Compared  to  the  curve  in  Fig.  3 
with  b/a  - 100,  the  present  curve  lies  slightly  lower. 


5. 2 Corona  Effect  on  Large  Signal  Propagation 

Consider  a semi-infinite  straight  wire  of  radius  a hung  over  a perfectly 

conducting  ground  plane  at  a height  h,  as  illustrated  in  Fig.  9.  It  is  taken 

to  stretch  from  * - 0 to  «.  Let  the  wire  be  quiescent  for  t < 0.  If, 

starting  at  t “ 0,  a voltage  difference  V (t)  is  applied  between  the  free  end 

8 

of  the  wire  and  the  ground  by  a generator,  the  signal  will  propagate  down  the 
line  to  infinity.  To  a good  approximation,  the  propagation  is  described  by  the 
pair  of  transmission  line  equations: 


.31  11  . r 

3z  " 3t  ’ 3z  3t 


(113) 


where  V is  the  shunt  voltage,  I the  series  current,  and  L and  C are, 
respectively,  the  inductance  and  capacitance  per  unit  length  of  line.  For 
h » a,  one  has 


o /.  /2h. 


C » 4irc  [4in(— )] 
o a 


-1 


(114) 


If  the  signal  exceeds  the  breakdown  voltage  for  the  wire,  coronas  will 

be  formed  along  its  path.  Current  will  flow  from  the  wire  to  the  ground,  forming 


a leakage  and  a loss.  In  this  event,  the  transmission  line  equations  (113) 
must  be  corrected  by  a corona  shunt  current  per  unit  length  1(V): 


(114) 


In  the  following  it  will  be  assumed  that  1(V)  can  be  given  its  d-c  value  from 
Eq.  (109).  This  is  a nonlinear  function  of  V.  With  the  introduction  of  a 
corona  shunt  conductance  per  unit  length  G such  that  i > GV,  the  equivalent 
circuit  for  Eqs.  (114)  is  represented  in  Fig.  10. 

Elimination  of  I from  (114)  yields  the  second-order  nonlinear  wave 
equation  for  V: 


2 " 2 3t. 

3z  C at 


(115) 


where 


(116) 


is  the  velocity  of  light  in  free  space.  For  i ■ 0,  as  when  V < V^,  Eq.  (115) 
reduces  to  the  standard  linear  wave  equation  whose  solution  is  of  the  form 


■ V(t  - z/c) . 


(117) 


This  represents  a voltage  disturbance  propagating  undeformed  and  undamped  with 
the  speed  of  light  in  the  positive  z-dlrectlon.  For  1^0,  this  simple  form 
is  no  longer  adequate.  Instead,  one  must  ass\mie 


V - V(z,  t-  z/c). 


(118) 


This  added  direct  dependence  on  z accounts  for  the  distortion  of  the  signal 
due  to  the  formation  of  coronas  along  the  line.  It  is  convenient  to  change  the 
independent  variables  from  the  set  (z,t)  to  the  set  (z,t)  with 
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Figure  10.  Equivalent  circuit  for  transmission  line,  with  corona  loss  represented 
by  nonlinear  shunt  conductance  per  unit  length  G(V). 


45 


i 

y 

i 


T ■ t - z/c. 


(119; 


whereby 


9 I ^ 9 I ^ 9 I 9 I ^ 9 I 

*t  92  c 9t  ‘z’  9t  ‘z  9t  'z‘ 


(i2o; 


Substituting  this  into  (115) , one  obtains 


9^V  2 9^V  ^ 

- 2 “ c 9t9z  9t 

oZ 


(i2i: 


In  the  following  analysis,  the  corona  loss  will  be  assumed  small.  Following 

14 

a method  of  Ilinova  and  Khokhlov  , one  neglects  the  second-order  derivative  with 
respect  to  z in  Eq.  (121).  This  is  because  if  the  corona  distortion  is  small, 
the  solution  is  approximately  of  the  form  (117)  and  the  direct  dependence 
on  z will  be  weak,  so  that  only  the  first-order  derivative  with  respect  to 
z need  be  retained.  Then  Eq.  (121)  becomes 


_L  rl  9V 

9t  ^c  9z 


+ Li) 


0, 


(122) 


with  the  solution 


C(z) 


(123; 


where  C is  function  of  z alone.  Eq.  (123)  holds  for  all  t.  If  the  line  is 

unexcited  for  t < 0,  one  must  conclude  that  C(z)  ■ 0.  Introducing  the 

characteristic  impedance  Z of  the  line,  such  that 

c 


Z - cL 
c 


c 

o 


‘na  , 


(124) 


one  can  rewrite  (123)  as 


} $r.t: 


(12 


3V  1 

^+2  =0- 


For  a given  functional  dependence  of  i on  V,  the  integration  of  this  equation 
is  straightforward.  One  assumes  here  for  simplicity  the  parabolic  form  (55): 


V < V 


icV(V-V  ) V > V 
c c 


The  constant  k can  be  adjusted  to  fit  Eq.  (109)  or  Fig.  8.  Then  the  solution 
of  (125)  is 


V(z,T) 


1 -g(T)exp(-^  Z <V  z) 
z c c 


V < V 


V > V 


where  f(T)  and  g(T)  are  some  functions  of  r alone,  to  be  determined  by 
boundary  conditions. 

Suppose  the  generator  in  Fig.  9 puts  out  a voltage  pulse  of  the  double- 
exponential  form: 


* vJe-“^-e‘®")9(t), 


and  suppose  that  is  sufficiently  large  that  ^g(^)  exceeds  the  breakdown 

voltage  for  t'  < t < t",  as  shown  in  Fig.  11.  The  unknown  functions 

f(T)  and  g(T)  in  (127)  will  be  determined  by  matching  the  solution  V(z,t) 
with  the  boundary  value  (128)  at  the  driving  point  z » 0 of  the  line.  At  this 
point,  T is  equal  to  t,  and  (127)  becomes 


Equating  this  with  (128)  yields  the  complete  solution  for  general  z 

and  t: 


0 < T < t' , 
t'  < T < t". 


t"  < t < *. 


(130) 


This  result  is  easy  to  interpret.  The  portions  of  the  pulse  in  Fig.  11  below 
t'  and  above  t",  being  sub-critical,  will  propagate  unaltered  and  unattenuated 
down  the  line  wi.th  the  speed  of  light.  The  portion  between  t'  and  t"  will 
be  distorted  by  corona,  and  will  be  leveled  down  to  at  a sufficiently  large 

distance  z from  the  voltage  input  for  which  the  exponential  factor  in  (130) 
is  small.  The  velocity  of  propagation  and  the  width  of  this  portion  are  unaffected. 
This  result  can  readily  be  extended  to  any  general  pulse  input  other  than  (128): 
all  supercritical  portions  of  the  pulse  will  eventually  be  chopped  down  by  corona 
to  the  critical  value. 

The  order  of  magnitude  of  the  propagation  distance,  at  which  corona 

attenuation  is  appreciable,  can  be  estimated  from  the  damping  constant  Z icV  in 

—6  ^ ^ 

(130).  Typically  this  quantity  has  a value  '>'10  /meter.  Thus  the  pulse  has 
to  travel  some  1000  kilometers  before  corona  attenuation  becomes  severe.  The 
reason  for  the  smallness  of  the  effect  Is  not  hard  to  find.  The  corona  onset 
voltage  is  usually  very  high.  At  this  voltage,  the  series  current  in  the  transmission 
line  is  proportionately  large,  being  normally  of  the  order  of  1000  A.  The  corona 
current,  on  the  other  hand,  is  typically  of  the  order  of  10  A/meter. 

Thus,  in  one  meter  of  line,  the  current  loss  is  only  about  one  part  in  a million. 

The  extreme  smallness  of  the  attenuation  amply  justifies  the  neglect  of  the 
second-order  partial  derivative  with  respect  to  z in  (121) . 


SECTION  6 


CONCLUSIONS 


This  report  comprises  a study  of  the  physics  and  the  mathematics  of  the 
electrical  corona  on  a long  wire.  The  fundamental  physical  processes  of  collision 
ionization,  ion  drift  and  diffusion,  which  underlie  much  of  the  gas  breakdown 
phenomena,  were  discussed  at  considerable  length.  They  were  transcribed  into 
precise  mathematical  language  In  terms  of  a set  of  coupled  nonlinear  partial 
differential  equations.  A number  of  boundary-value  problems  for  these  basic 
equations  were  worked  out.  The  calculations  for  the  d-c  and  microwave  breakdown 
fields,  and  for  the  d-c  voltage-current  characteristic,  were  found  to  be  in 
good  agreement  with  experiment.  These  successes  provided  solid  support  to  the 
correctness  of  the  basic  equations. 

With  the  exception  of  the  transmission-line  problem  in  Section  5,  the 
boundary-value  problems  solved  in  this  study  are  time-independent  problems. 

This  is  not  to  suggest  that  the  time-dependent  problems  are  not  of  primary 
interest,  but  rather  that  they  are  by  far  harder  to  tackle.  The  mathematical 
equivalent  of  a transition  from  a time-independent  corona  phenomenon  to  a 
time-dependent  one  is  a passage  from  nonlinear  ordinary  differential  equations 
to  nonlinear  partial  differential  equations.  The  range  of  mathematical  techniques 
available  today  for  handling  the  latter  type  of  equations  is  severely  limited. 

And  yet  the  study  of  many  important  aspects  of  the  electrical  corona  will  involve 
time-dependent  problems. 

In  Subsection  5.2,  the  d-c  corona  voltage-current  characteristic  was 
applied  to  calculate  the  corona  effect  on  pulse  propagation  along  a transmission 
line.  It  is  to  be  noted  that  although  the  problem  is  time-dependent,  the  treatment 
of  the  corona  itself  is  quasi-static.  This  procedure  is  reasonable  only  if  the 
shunt  voltage  at  a fixed  point  on  the  line  changes  but  little  during  the  time  it 
takes  for  the  corona  current  to  establish  a steady  flow.  The  latter  time  is 
at  least  as  long  as  the  transit  time  for  an  ion  to  travel  from  one  conductor 
to  another.  An  ion  in  air  typically  picks  up  a velocity  of  1 cm. /sec.  in  a 
field  of  1 volt/cm.  At  breakdown,  a representative  ion  velocity  is  about 
1 km. /sec.  If  the  conductor  separation  is  1 meter,  the  transit  time  is  about 

-3 

10  second.  One  can  surmise  that  if  the  width  of  a pulse,  or  the  period  of 
a continuous  wave,  is  smaller  than  this  value,  the  d-c  result  will  be  inapplicable. 


A time-dependent  treatment  of  the  corona  current  must  be  considered. 

Another  problem  in  which  time-dependence  Is  essential  Is  that  of  the 
corona  effect  on  the  traillng-wire  antenna.  During  operation,  the  greater  part 
of  this  antenna  Is  hung  In  mld-alr,  far  away  from  other  conductors.  There  is 
no  d-c  equivalent  of  this  configuration,  since  the  d-c  corona  current  always 
flows  between  two  electrodes.  The  fact  is  that,  under  a continuous-wave  excitation, 
the  antenna  is  alternately  anode  and  cathode.  The  ions  in  the  corona  oscillate  in 
and  out  radially  near  the  wire  surface,  more  or  less  in  phase  with  the  wave. 

If  there  is  a constant  corona  ohmic  loss  per  cycle,  then  the  power  loss  due  to 
corona  is  directly  proportional  to  the  wave  frequency,  as  was  experimentally 
found  by  Peek  in  Eq.  (3).  At  high  frequencies,  the  loss  is  potentially  very 
large.  Obviously  a theoretical  attempt  to  deduce  Peek's  empirical  result  entails 
the  solution  of  a time-dependent  problem. 

The  analysis  of  the  corona  effect  on  aircraft  antennas  is  a challenging 
theoretical  problem.  Its  difficulty  Is  enhanced  by  the  scarcity  of  antenna  data 
to  guide  the  calculations.  This  preliminary  study  has  been  directed  at  identifying 
the  basic  physical  processes  and  the  basic  mathematical  problems  Involved.  It 
has  been  concluded  that  Eqs.  (9)  through  (13)  are  the  basic  Irreducible  equations 
for  the  corona  phenomenon.  Their  essential  correctness  has  been  established 
by  the  good  agreement  between  calculation  and  measurement  in  the  quasi-static 
limit.  The  next  stage  in  the  overall  corona  effort  will  Involve  the  study  of 
time-dependent  solutions  of  the  basic  equations.  These  will  describe  both  the 
EMP  and  CW  excitations  of  long-wire  antennas.  The  chief  quantities  to  come  out 
of  the  calculations  will  be  the  corona  charge  and  current  densities  around  the 
antennas  as  a function  of  space  and  time,  and  of  the  sources  of  excitation. 

From  these  results  one  can  analyze  such  corona  effects  as  antenna  signal  distor- 
tion, power  loss,  and  the  corona  modification  of  the  coupling  between  aircraft 
and  antennas. 
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